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Abstract. We performed optical-pump terahertz-probe measurements of a Mott 
insulator YTiC>3 and a band semiconductor Si using a laser diode (1.47 eV) and a 
femtosecond pulse laser (1.55 eV). Both samples possess long energy-relaxation times 
(1.5 ms for YTiC>3 and 15 /is for Si); therefore, it is possible to extract terahertz 
complex conductivities of photoinduced carriers under equilibrium. We observed 
highly contrasting behavior - Drude conductivity in Si and localized conductivity 
possibly obeying the Jonscher law in YTiC>3. The carrier number at the highest 
carrier-concentration layer in YTiC>3 is estimated to be 0.015 per Ti site. Anisotropic 
conductivity of YTiC>3 is determined. Our study indicates that localized carriers might 
play an important role in the incipient formation of photoinduced metallic phases in 
Mott insulators. In addition, this study shows that the transfer-matrix method is 
effective for extracting an optical constant of a sample with a spatially inhomogeneous 
carrier distribution. 
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1. Introduction 

Recent discoveries of photoinduced metallic phases in several Mott insulators [H El [31 
HJ [5l [6] made us consider the strongly correlated electron physics from a new point 
of view. However, the optical properties of photoinduced carriers in Mott insulators 
are not well understood even at low excitation densities. Understanding these optical 
properties is a prerequisite to understanding the incipient creation of metallic phases. 
Drude response by itinerant carriers is observed in the case of band semiconductors with 
low excited-carrier density (10 14 ~10 16 cm -3 ) [3, El El ED]- The comparison between the 
optical properties of photoinduced carriers at low excitation densities in Mott insulators 
and those in band semiconductors would be important in gaining deeper insight into 
strongly correlated electron physics. 

The detailed nature of various carrier conductions, exhibiting Drude or hopping 
conduction, can be well characterized in the terahertz (THz) regime [HI [121 [El EH]- 
THz time-domain spectroscopy (THz-TDS) is a powerful tool for analysing terahertz 
conductivity a(ou) (= o~i{uj) + ^(cu)). The remarkable advantage of THz-TDS is its 
simultaneous determination of both the real and imaginary parts of without using 
the Kramers-Kronig transformation [13]. The coherent nature of the THz pulse is also 
utilized to investigate photoinduced by, for instance, optical-pump THz-probe 

(OPTP) studies. Recent progress in THz technologies [T6l [T7] and the methodology 
of analysis [181 [T9l [20] in OPTP experiments enable evaluation of transient (t(u) in 
many substances, such as semiconductors, high-T c superconductors, liquids and organic 

materials^ 12111221123121123120112312^ 

Because the OPTP method essentially detects non-equilibrium processes, such as 
surface recombination and carrier diffusion, the time dependence of a complicated spatial 
carrier distribution must be considered [8l [32]. This difficulty is avoided by using thin- 
film samples [2T1 [221 [231 EH EHl [31] where the optical pump pulse penetrates, and by 
analysing the photoinduced phase as a layer with a homogeneous cr{ui) [H E21 EH 1281133]. 
Furthermore the extraction of varying quickly compared with the pulse width 
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of THz probe pulse, only seems possible within some restricted conditions [TS1 ESI |20j . 
Therefore, to analyse photoinduced for a wide variety of Mott insulators, free 

from the restrictions in sample preparation and analysis, we initially examined a nearly 
equilibrium state of photoexcited bulk material with a longer energy-relaxation time, 
t. The photoexcitation of a material with longer r creates a quasi-equilibrium state 
averaging over various non-equilibrium processes, making it easy to obtain the optical 
constants of the highest carrier-concentration region in the material. 

Our analysis also required extracting a(u) of materials with inhomogeneous carrier 
distributions in a more rigorous manner. One good candidate for accomplishing this 
is the transfer-matrix method, which expresses inhomogeneous carrier distribution by 
a multi-layer system. This method has been briefly commented on in the literature [8]. 
The promise that the transfer-matrix method can incorporate inhomogeneity is seen 
in the analyses of reflectivity in optical pump-probe studies[3U [35]. Although it is 
a versatile method, its effectiveness has not been thoroughly discussed, especially in 
THz-TDS studies. 

In this study, we found that r of a Mott insulator YTi03 with Mott gap of 
approximately 1 eV[36] is 1.5 ms at 1.47 eV photoexcitation. We characterized the 
photoinduced &(u) by comparing it with that of a band semiconductor Si with a bandgap 
of 1.1 eV[37]. We also present a more detailed discussion of the transfer-matrix method. 

2. Experimental method 

Single-crystalline samples of YTiOs, with the orthorhombic perovskite GdFe03-type 
structure, were grown by the floating zone method [38]. The Si sample was commercial 
high-resistivity Si. 

Photoconductivity was measured to assess r. The light emitted from a multimode 
continuous wave (CW) laser diode (LD) with photon energy of 1.47 eV was modulated 
and used for illuminating the sample under an electric field of about 0.3 kV/cm. The 
photocurrent I p h flowing through a 100 Q resistance, connected in series with the sample, 
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was lock-in detected. 

OPTP experiment was performed by a transmission THz-TDS system described in 
detail elsewhere [39, 40J. The thicknesses of platelet samples are 420 /mi for YTi03 and 
512 fim for Si. The optical pulses were generated by a mode-locked Ti-sapphire laser with 
a repetition rate of 76 MHz and central wavelength of 800 nm. Both the THz emitter 
and detector were low-temperature grown GaAs photoconductive antennas. Si lenses 
were attached to the antennas to enhance the emission power and collection efficiency 
of THz pulses. The THz spectral range in this experiment was between 0.5 and 8 meV. 
The THz-wave-emission sides of the samples were photoexcited by the multi-mode CW 
LD with an incident angle of 45°. The pump-beam power was 0.8 W for Si and 1.2 W 
for YTiCV For YTiC>3, the polarization of THz electric field Ethz is parallel to the 
6-axis. The beam diameter of the LD light was about 8 mm and larger than that of the 
THz probe pulse, which is energy dependent (e.g. 2 mm at 2 meV and 1 mm at 4 meV). 
The fluence rate was 1.6 W/cm 2 for Si and 2.4 W/cm 2 for YTiC>3, respectively. The 
temperature rise resulting from the thermalisation by photoexcitationjjl] is estimated 
not to exceed 1 K. 

OPTP measurements were also performed with the 1.55 eV optical-pump pulses 
split from the Ti-sapphire laser to investigate the nature of conduction of carriers induced 
by light with photon energy larger than that of the CW LD. This experiment studied 
the anisotropy of photoinduced &(uj) in YTiOs. The optical-pump pulse power was 230 
mW. The optical-pump pulse beam diameter was about 2 mm (fluence: 96 nJ/cm 2 ) - 
slightly smaller than that of the THz probe pulse below 2 meV. The Ethz was applied 
along either the b- or c-axis, maintaining the polarization of the optical-pump pulse 
electric field Ei 55e y parallel to the b- or c-axis. 

All measurements were performed at room temperature. 
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3. Results and Discussion 

Figure 1 shows the modulation-frequency dependence of I p h- Increasing the modulation 
frequency causes I p h to decrease, according to ( 1+ (f r )2- ) i where B is the proportional 
coefficient. The obtained B and r are listed in the figure. Longer r (~ 15 fis for Si and 
~ 1.5 ms for YTi0 3 ) indicates that the samples are in quasi-equilibrium photoinduced 
state during the THz-TDS measurements, r of YTi0 3 is almost excitation-intensity 
independent; this implies that a thermal effect does not dominate the relaxation process. 
Although the direction of I ph and of electric field of the CW LD pump are not identified 
accurately in YTi0 3 , a huge anisotropic r that depends on the direction of I ph and the 
polarization of the excitation light is not anticipated. The OPTP results shown later 
support this (see Fig. 7). 

The temporal evolution of E THz transmitted through Si and YTi0 3 are shown 
in Fig. 2(a) and 2(b), respectively, with and without LD excitations (1.47 eV). 
Photoexcitation attenuates both THz waves implying THz wave absorption is by the 
photoinduced carriers. THz energy dependence of transmission T{uS) and phase shift 
A(f)(uj) are obtained by Fourier transformation of THz waves, as shown in Fig. 3(a) 
and 3(b). They are calculated with the equations T(tu) = and A0(cj) = 

4>{ui) — (f) re f(ui), where E(uj){E re f{uj)) and <p(uj){(j) re f{uj)) are the Fourier transformed 
amplitude and phase, with and without excitation, respectively. THz wave absorption 
by photoinduced carriers is responsible for T{uj) decreasing below 1 for both samples, 
but the two exhibit different energy dependencies. As the THz energy increases, T(uo) of 
Si approaches 1, while that of YTi0 3 gradually decreases. The opposite sign of A0(u;) 
of the two samples strongly indicates different fundamental conduction mechanisms of 
the photoinduced carriers. Negative (positive) A0 roughly means the refractive index 
is reduced (increased), compared with an unexcited state, which influences the negative 
(positive) real part of the dielectric constant of photoinduced carriers. As explained 
below, these results indicate a metallic nature below a plasma frequency in Si, and a 
localized nature, such as hopping carriers, in YTi0 3 . 



Terahertz conductivity of localized photoinduced carriers in Mott insulator YTiO^ 7 
Before showing the photoinduced (j(uj) of Si and YTi0 3 , we mention the detailed 
procedure of the transfer-matrix method. A spatially inhomogeneous distribution of 
photoinduced carriers is initially regarded as exponentially decaying. Subsequently a 
non-exponentially decaying distribution is introduced. The photoinduced phase with 
an exponentially decaying carrier distribution is divided into many thin slabs (see Fig. 
4). Each slab is supposed to have a uniform complex refractive index fij whose value is 
set to reproduce the exponential decay of the photoinduced carrier concentration. The 
transfer-matrix of each slab is described by 



M, 



1 cos 5 i sin 5/rij ^ 
i ihj sin 5 cos 5 j 



(1) 



where j is the number index of the slab, 5 = ^r-hjd, A is the incident THz wavelength 
in vacuum, and d is the slab thickness |42[ |4"3] . The fij of each slab is calculated from 
the complex dielectric constant e using, 

n] = Ino + e sur exp > ( 2 ) 

where e no is e without excitation, i sur , the parameter to be optimized in this analysis, 
is e resulting from carriers at the photoexcited surface of the sample, Zj (= j x d) is the 
depth from the photoexcited surface into the sample along the THz wave propagation, 
and d p is the optical penetration depth. For Si, frequency independent e no [H] of 11.7 is 
used. For YTiC>3, e no is experimentally determined by the THz-TDS measurement and 
is weakly energy dependent (e.g. 16.5+0.4i at 2 meV and 17+0. 4i at 4 meV). d p of Si at 
1.47 eV is determined to be 8.4 /xm using the absorption coefficient from an optical data 
handbook[45j. That of YTiC>3 at 1.47 eV was calculated to be 0.22 //m from the reported 
reflectivity spectra[46] (0.05-40 eV) combined with the Kramers-Kronig transformation. 
Then the total matrix M t is described as 

M t =f[M v (3) 

j=k 

where k is the total number of photoexcited slabs. Finally, the THz complex 
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transmission is given by 

rp/ \ C\i( \\ tw/LD (Q + iP)wo/LD 

T(uj) exp(zA0(w)) = —!■ = — - — , (4) 

two/LD (,W + 1 -l)w/LD 

where t w /LE> and t^o/LD mean the THz wave transmission with and without excitation, 
respectively, Q = Re((M tll +M tl2 )^ +M m +M t22 ) and P = Im((M m +M tl2 )^ + 
M t2 i+M t22 ). The fitting of experimental T(u) and A<p(u) following the above mentioned 
procedure provides cr(o>) resulting from carriers at the photoexcited surface through 

a(u) = iue e sur (uj), (5) 

where eo is the vacuum permittivity. Note that the convergence of transmission is 
checked carefully by decreasing the thickness or by increasing the number of slabs. A 
thickness of the photoexcited phase (= k x d, fc=100) 5 to 10 times thicker than d p is 
typically employed. 

Figure 5 shows photoinduced of Si and YTi03. cr{uj) of Si can be interpreted 
by the Drude model as 

V = — ^ ( 7 ) 
where n c is the carrier density, \i is the mobility, T is the carrier collision rate and m* is 
the effective mass. The photoexcitation introduces both electrons and holes; therefore, 
the tentatively assigned m* value is 0.26mo for electrons and 0.37m for holes[H], where 
mo is the free-electron mass. Hereafter, fi of each carrier is denoted as /i e for electrons 
and fih for holes. We have considered the following two cases, neither of which can be 
excluded at the present stage. One is the two-carrier model of electrons and holes. The 
other takes only electrons into consideration, assuming that holes with heavy m* do not 
contribute to a{oj). The solid lines in Fig. 5(a) represent the calculated a(u) for the 
two-carrier model and the broken lines represent the electron-only model. The curves 
are in agreement with the experimental a(uj). This suggests that the itinerant carriers 
are certainly photogenerated in Si. The obtained /x e and /x^ are 2410 (±210) cm 2 /Vs 
and 500 (±90) cm 2 /Vs for the two-carrier model, and \x e is 1820 (±100) cm 2 /Vs for the 
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other model. They are roughly consistent with the literature values [57], but it is to be 
noted that \x e in both models might be larger than the predicted ones. The ambiguity 
of m* may be responsible for this deviation. 

The most striking feature in a{uj) of YTiC>3 is the negative o~i. It suggests an 
existence of localized carriers [27 | |39 | |4"T] . which is very different from Si. The localization 
may arise from the on-site strong Coulomb interaction between 3d electrons in YTiC>3. 
To explain a(u), we used the empirical Jonscher law [IS], which expresses cr(a>) of many 
materials with hopping carriers. The Jonscher law is given by[4"8l H9] 

<ti(w) = (J dc + Au s , (8) 
a 2 (ou) — -Au s tan— , (9) 

where is the DC conductivity, A the proportional coefficient and s is restricted 
between and 1. As shown in Fig. 5(b), the solid curves from the Jonscher law 
seem to agree with experimental a(u). In the solid curves, s, a^c and A are 0.95, 
235 (±10) n^cm- 1 and 2.40 (±0.05) x lO" 11 fi^an^s - 95 , respectively. The allowed 
s ranges from 0.91 to 0.99, and corresponding o^c and A are 210 (±10) r2~ 1 cm _1 and 
1.41 (±0.03) xlO" 10 fi^cm-V- 91 , and 260 (±10) fi^cm" 1 and 1.43 (±0.03) x 10" 12 
r2 _1 cm _1 s ' 99 , respectively. 

Note that o(uS) can be also fitted by a two-component model, such as the Drude- 
Lorentz. The estimated photoinduced carrier number at the surface layer is about 0.015 
per Ti site. Photoexcited YTiOa with the derived carrier density would be equivalent 
to chemically hole-doped Y 1 _ x .Ca :r Ti03 with x much less than 0.1 given in Ref. 36. 
The o~x(ou) spectrum of Yi^Ca^TiOs in this composition region is very different than a 
Drude response. Therefore, it would be difficult to expect a Drude component to exist. 
Clarifying this point might require broadband spectroscopic information obtained under 
photoexcitation or the temperature dependence of a{uj). 

Since long relaxation times, r, are observed in both samples, a diffusion or a 
surface-recombination process, making the carrier distribution a non-exponential decay 
type, must be considered, and the analysis method modified. The carrier number n(z) 
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along the THz wave propagation in a quasi-equilibrium state is obtained using a one- 
dimensional diffusion equation [50] as follows: 



where n(z, t) depends on the time t and the position z along the THz wave propagation, 
and 5(t) is the 5-function. D is the diffusion coefficient and given by 

D = ^El i (11) 
\e\ 

where \j \iu is equal to l//x e + 1/ fih, &b is the Boltzmann constant and T s is the sample 
temperature equal to 300 K. The solution[50] of eq. (ITUI) is 

z 2 \ I I 



n(z, t) = exp 



ADt / 2 



f v -v s J \ \D 2VDi) I V r 
where v s is the surface recombination velocity and f(z) is related to the error function 
by f(z) = exp(z 2 )(l— erf(z)). The carrier number in the quasi-equilibrium state requires 
the integration of n(z, t) with respect to t, 

/•OO 

n(z) = / n(z,t)dt. (13) 



o 



Therefore, with the assumption of a conduction model and the knowledge of n(z) 
determined by appropriate /i^ and v s , T(uj)e lA ^^ can be calculated using eq. (@J. 
In this case, eq. ([2]) is replaced by 



where n max and i max are n(z) and e of the highest carrier-concentration layer, 
respectively. 

After the determination of \iu of Si (= 346 cm 2 /Vs) using the literature values [37]. 
v s is varied between lxlO 4 cm/s and lxlO 6 cm/s. Representative n(z) normalized at 
n m ax are shown in the inset of Fig. 6(a). The n max is observed around 1~2 /xm. Assuming 
that both electrons and holes obeying the Drude conductivity are responsible for a(u), 
T(uj)e lA ^ u '' > is confirmed as being consistent with experimental data for both n(z) (see 
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Fig. 6(a)). The estimated n c is 5.2 (±0.3) xlO 16 cm -3 and is comparable to that obtained 
by the previous model. This indicates that, at the highest carrier-concentration layer 
almost the same n c can be obtained, irrespective of the carrier distribution decay type. 
For YTiC>3, both v s and \iu are unknown parameters. The wide range sweep of v s and \iu 
gives various n(z) curves as depicted in the inset of Fig. 6(b) with peak positions around 
0.1 /im. For each n(z), the experimental T(cu)e lA ^^ is well reproduced by the Jonscher 
law, where s is restricted within the same range obtained in Fig. 5(b) (0.91~0.99). 
Typical examples are shown in Fig. 6(b) with s of 0.95. The other parameters (a^c 
in fi _1 cm _1 and A in f2 _1 cm~ 1 s ' 95 ) for the dotted-solid, solid and broken lines are 
125(±10) and 1.40(±0.05)x 10 -11 , 125(±10) and 1.2(±0.1) xlO" 11 , and 90(±10) and 
1.0(±0.1)xl0 -11 , respectively. o~\(uj) and o~i{oS) calculated from the parameters are half 
to two-thirds of those in Fig. 5(b). Thus, for YTi03, the a(a>) extracted from the 
model with exponentially-decaying carrier distribution roughly represents the highest 
carrier-concentration layer in the model using eq. ([141) . 

The photoinduced carrier number at the highest carrier-concentration layer in 
YTiC>3 is calculated as 0.015 per Ti site[51]. It can be proposed, therefore, that 
a phase with localized carriers would emerge initially at the photogeneration of the 
metallic phase in Mott insulators. Photoexcitation creates both electrons and holes, 
which differs from chemical doping, and a comparison of &{uS) between photoexcited 
YTiC>3 and a hole-doped Yi-^Ca^TiOs is discussed. The absolute value of cri(u) of 
photoexcited state might be much larger than that of corresponding Yi-^Ca^TiOs if 
the extrapolation of o~i(u) in Y 1 _ a .Ca ;c Ti03 is carried out toward the THz energy. The 
preservation of spectral weight implies that the localization energy of photoinduced 
carriers would be much lower than for holes in Y^^-Ca^TiC^, even if both holes and 
electrons contribute to cr(u;) in photoexcited YTiC>3. As it is not clear that the large 
difference in localization energy originates from only holes in such a low-carrier system, 
it is plausible that electrons with small localization energies also contribute to a(u). 
Therefore of photoexcited YTiC>3 would be supported by bound electrons as well 
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as holes. 

In a halogen-bridged Ni one- dimensional chain compound [Ni(chxn)2Br]Br2 (chxn = 
cyclohexanediamine), which is compared with YTi03 composed of a three-dimensional 
Ti network, the localized o\ is determined at a lower excitation density^]. Despite being 
in a different energy region, carrier localization in photoexcited Mott insulators at low 
excitation densities may be the general phenomenon, irrespective of the dimensionality. 

The fact that one-dimensional Mott insulators, such as [Ni(chxn)2Br]Br2[3] and 
Sr 2 Cu0 3 [52], exhibit r in the order of pico-seconds may suggest that dimensionality is 
a decisive factor of r. 

Figure 7 shows T(u) and A<f>(u) obtained by OPTP experiments using a 
femtosecond-pulse laser (1.55 eV) for YTi03. Since the period of optical-pump arrival 
time (13 ns) is much shorter than r, the photoinduced carriers are also in quasi- 
equilibrium state. In both polarizations of Ei 55e y, it is found that the degree of variation 
from the unexcited state in THz wave amplitude and phase is larger for Ej<H Z \\b within 
the measured THz energy range. This implies that the absolute values of o~\ and o~i 
for EthzWo are larger than those for Ej<hz\\ c - The anisotropy would reflect the crystal 
symmetry of YTiC>3 or the 3<i-orbital state at Ti site. Comparing with the OPTP 
results to those using the CW LD, the a(u) does not seem to depend strongly on the 
optical-photon energy. 

4. Summary 

We have optically characterized photoinduced carriers of Mott insulator YTiC>3 at low 
excitation densities in the THz regime by OPTP measurements, and compared the 
experimental results with those for band semiconductor Si. The r of the photoinduced 
carriers in YTi0 3 is about 1.5 ms. The inhomogeneous carrier distribution along the 
THz wave propagation can be treated accurately using the transfer-matrix method. This 
method successfully determined of the highest carrier-concentration layer under the 
quasi-equilibrium states. YTi03 shows localized a(uj), possibly with the Jonscher law, 
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whereas Si exhibits the Drude response. Anisotropic 5{uj) in YTi0 3 is determined. Our 
study demonstrates that localized carriers might play an important role in the incipient 
formation of metallic phases in photoexcited Mott insulators. Although the exact origin 
of the localization in YTi0 3 remains an open question, THz-TDS under photoexcitation 
with another photon energy or for another Mott insulator might provide the answer. 
We note here that a preliminary THz-TDS experiment of YTi0 3 excited by a CW LD 
of 1.9 eV also leads to localized o(uo). 
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Figure 1. Photocurrent as a function of modulation frequency for Si and YTiC>3. The 
solid lines are B/(l + (cut) 2 ) calculated with B and r denoted in the figure. 
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Figure 2. Time evolution of Ethz transmitted through (a) Si and (b) YTi03 with 
(solid lines) and without (broken lines) excitations (1.47 eV). The polarization of Ethz 
is parallel to the 6-axis in YTiC>3. The excitation power is 0.8 W for Si and 1.2 W for 
YTi0 3 . 



Terahertz conductivity of localized photoinduced carriers in Mott insulator YTiOs 19 



(a) Si 
1.5 L 



t 1 1 r i" 



***** 



E 1 



0.5 



'*»** 



so 
o 

-50 



-100 -£ 
-150 <| 



200 



(b) YTiO, 



1,5 



E 1 



****** 



-250 



40 



CD 
T3 



0.5 - 



_t l I L 1_ 



-40 < 



-so 



I 2 3 4 5 G 7 S 
Energy (meV) 



Fig.3 



Figure 3. THz energy dependence of transmission and phase shift for (a) Si and (b) 
YTiC>3. The polarization of Ethz is parallel to the 6-axis in YTiC>3. 
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Figure 4. Illustration of a photoexcited sample to explain the transfer-matrix method. 
The coloured area demonstrates an exponentially-decaying carrier density. The area 
is divided into many thin slabs with thickness d. The total number of thin slabs is k, 
and the transfer-matrix of each slab is expressed as Mj (j = 0, 1, k). The depth 
from the excited surface into the sample along the THz wave propagation is denoted 

clS Zj . 
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Figure 5. THz complex conductivity of (a) Si and (b) YTi03 at the photoexcited 
surface. The real and the imaginary part of conductivity correspond to u\ and cr 2 , 
respectively. The solid and broken curves in (a) are <j(cu) calculated by the Drudc 
models with n c =6.6 (±0.3) xlO 16 cur 3 , /i e =2410 (±210) cm 2 /Vs and /x h =500(±90) 
cm 2 /Vs for the solid curves (two-carrier model of electrons and holes), and with 
n c =9.3(±0.4)xl0 16 cm~ 3 and /x e =1820 (±100) cm 2 /Vs for the broken ones (only 
electrons under consideration), respectively. In the two-carrier model, the same n c 
is assumed for each carrier. The solid curves in (b) are calculated <j(ui) using the 
Jonscher law with a dc =235 (±10) fi^cm -1 , A =2.40 (±0.05)x lO" 11 ^cm^s - 95 
and s =0.95. 
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Figure 6. THz energy dependence of transmission and phase shift for (a) Si and 
(b)YTi03 analysed by the model with non-exponential carrier-distribution decay. The 
evaluated parameters v s and \±u are listed in the figure. In (a), Drude conductivity 
characterized by n c —5.2 (±0.3) xlO 16 cm~ 3 and the literature values of \i e =1500 
cm 2 /Vs and \ih =450 cm 2 /Vs is assumed. The conduction model employed in (b) is 
the Jonscher law (s=0.95), where Ode in il _1 cm _1 and A in Q _1 cm _1 s 095 are 125 
(±10) and 1.40 (±0.05)xl0- n , 125 (±10) and 1.2 (±0.1)xl0- n , and 90 (±10) and 
1.0 (±0.1) x 10~ n for the dotted-solid, solid and broken curves, respectively. The insets 
of (a) and (b) show z dependencies of the carrier numbers normalized at n max . 
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Figure 7. THz energy dependence of transmission and phase shift of YTi03 obtained 
by OPTP experiments, using the femtosecond-pulse laser (1.55 eV) under £r.ffz||b and 
Ethz\\c for (a)Si.55 e v||b and (b)-Ei.55 e y||c The power of the optical-pump pulse is 
230 mW. 



